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Stability Analysis of Conformable 
Fractional Systems 
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Abstract 


In this paper, we investigate stability analysis of fractional differential 
systems equipped with the conformable fractional derivatives. Some stabil- 
ity conditions of fractional differential systems are proposed by applying the 
fractional exponential function and the fractional Laplace transform. More- 
over, we check the stability of conformable fractional Lotka-Volterra system 
with the multi-step homotopy perturbation method to demonstrate the effi- 
ciency and effectiveness of the proposed procedure. 


Keywords: Stability analysis; Asymptotical stability; Conformable frac- 
tional derivative; Lotka-Volterra system. 


1 Introduction 


Fractional differential equations are generalizations of classical differential 
equations of integer order that have recently proved to be valuable tools 
for the modelling of many physical phenomena, and have been the focus of 
many studies due to their frequent appearances in various applications, such 
as physics, biology, finance and fractional dynamics, engineering, signal pro- 
cessing and control theory [8,12]. Finding solutions to fractional differential 
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systems is rather complicated, consequently, the stability results of the frac- 
tional differential systems have been the main goal of the previous studies. 
For example, Matignon considers the stability of fractional differential sys- 
tems in control processing [11], Deng has studied the stability of fractional 
differential system with multiple time delays [6] and we have investigated the 
stability of fractional differential systems with Hilfer derivatives [15]. More 
novelty, authors in [3, 4,14] studied stability analysis of distributed order 
fractional differential equations with respect to the nonnegative density func- 
tion. There are many definitions of fractional derivatives, such as Riemann- 
Liouville, Grunwald-Letnikov and Caputo’s fractional derivatives [12], which 
these fractional derivatives do not satisfy the product rule, the quotient rule 
and the chain rule. In 2014, to overcome these and other difficulties, Khalil 
et al. [7] introduced a new simple well-behaved definition of the fractional 
derivative called conformable fractional derivative. This fractional deriva- 
tive is theoretically very easier to handle and also obeys some conventional 
properties that cannot be satisfied by the existing fractional derivatives, for 
instance, the chain rule [1]. However this fractional derivative has a weak- 
ness, which is the fractional derivative of any differentiable function at point 
zero. 


In this paper, we shall describe stability conditions for conformable frac- 
tional differential systems. In particular our analysis covers the linear con- 
formable fractional differential systems with commensurate order and in- 
commensurate order and the nonlinear conformable fractional differential 
system. At first, stability conditions will be established using fractional 
exponential function for linear conformable fractional differential systems, 
corresponding to this result we will derive asymptotic stability for the non- 
linear conformable fractional system. Then, we produce sufficient conditions 
for asymptotical stability of in-commensurate linear conformable fractional 
differential system by of the fractional Laplace transform and the fractional 
final value theorem. Finally, we present the multi-step homotopy perturba- 
tion method for obtain approximate analytical solution and stability of the 
conformable fractional Lotka-Volterra system to illustrate the validity of the 
results. 


2 Conformable fractional derivative 


Here, some basic definitions and properties of the conformable fractional cal- 
culus theory which can be found in [1,2,7] are presented. 


Definition 1. Let f : (0,0oo) > R, then, the conformable fractional deriva- 
tive of f of order a is defined as [7] 
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@G)@ sim Tt I) (1) 


e—+0 é. 


for all t > 0, a € (0,1). 
If f is a-differentiable in some (0,a), a > 0, and lim, tT (f) (t) exists, then 
t0 


by definition 


Ta (Ff) (0) = Jim, 1Ta (f) (. (2) 


The new definition satisfies the properties which are given in the following 
theorem. 


Theorem 1. Let a € (0,1], and f, g be a-differentiable at point t, then [7] 
(t) +Ta(af + 6g) =a1Ta(f)+61Ta(g), for alla,b eR. 
(tt) ¢Ta.(t") = pp th, for all we R. 
(itt) eTa(f9) = f eTo(g) +9 +To(f). 


(iv) +To(4) = geTa(D—f ital) In addition, if f is differentiable, then 


Tal f)(t) = te. 
In [1] T. Abdeljawad established the chain rule for conformable fractional 
derivatives as following theorem. 


Theorem 2. Let f : (0,co) — R be a function such that f is differen- 
tiable and also a-differentiable. Let g be a function defined in the range of f 
and also differentiable; then, one has the following rule 


+Ta(fog)(t) = (1 Taf) (9(t)) (e Tag) (t) 9(H)°™. (3) 


If t = 0, then 


1Ta(fog)(0) = lim, (: Taf) (g(6)) (eTag) (8) 9(6)°?. 


The fractional exponential function plays a very important role in the con- 
formable fractional differential equations. The fractional exponential function 
1lya@ 
e«' , where 0 < a < 1, is defined by the following series representation, 


OO pak 


akkl 
k=0 


1 
eat — 


Now, we list here the fractional derivatives of certain functions [7] 


(‘) Tales") = ea", 
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(it) ¢T(sin 4t*) = cos 44%, 
(itt) tT (cos 4t*) = — sin 40°, 
eer st. 


On letting a = 1 in these derivatives, we get the corresponding ordinary 
derivatives. 


Definition 2. (Fractional Integral) [7] Let a > 0 and t > a. Also, let f 
be a function defined on (a,t] and a € (0,1]. Then, the a-fractional integral 
of f is defined by, 


t 
x 
uas= [ ae, (4) 
if the Riemann improper integral exists. 


Theorem 3. (Integration by parts) [1], Let f,g : [0,b] > R be two func- 
tions such that fg is differentiable. Then 


b b 
| f(t) «To (9) (#) da(t) = fol’ | g(t) iTa(f)(t)do(t), (5) 


where da(t) = t?~'dt. 
It is interesting to observe that the a-fractional derivative and the a- 


fractional integral are inverse of each other as given in [7]. 


Theorem 4. (Inverse property). Let a € (0,1) and f be a continuous 
function such that ,Ig' f exists. Then 


Talele f(t) = f(t), for t>0. 


Definition 3. (fractional Laplace transform) [1] Let 0 < a < 1 and 
f : [0,co) > R be real valued function. Then the fractional Laplace trans- 
form of order @ starting from a of f is defined by, 


Ln{ H(t} = Fal) = f ” ens F(t) dat), (6) 
where da(t) = t*~1dt. 


Fractional Laplace transform for certain functions are presented as fol- 
lows [1] 
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oiegfitat, 38> 0. 
QP/@ 


© La{t?} = SryaF(1+4 4), s>0. 


Lofe=}= 7 s>l. 


2 s—l? 
> La{sinit*}=a, s>l. 
© Le{cosit*} = ah, s>1. 


Furthermore, using the properties of the fractional exponential function 
and integration by parts, we have 


LottTa (f) (t)} = 8 Fa(s) — f(0). (7) 


Next we prove a fractional version of final value theorem which will be 
useful is studying stability of conformable fractional systems. 


Theorem 5. (final value theorem) Let Fy(s) be the fractional Laplace trans- 
form of the function f(t). If all poles of sF.(s) are in the open left-half plane, 
then, 
lim f(t) = lim sFQ(s). (8) 
s—0 


too 
Proof. By using formula (7) we arrive at the following relation 
lim Lae (f) (O)} = lim (s Fa(s) ~ £(0)) = lim (s Fa(s)) — £0). (9) 


s—0 


Consider now the first term only. Since s is independent of t, the order of 
integrating and taking the limit can be interchanged 


lims_0 La{:Ta (f) (t)} = lims 40 fo" e7* = Ta (f) (t) da(t) 
= f° [limssoe**"] 1Ta(f)(#) date) (10) 
= fo° 1Ta (f) @) dott). 


Taking the integration in the last term of (10), implies that 


lim (f(t) — f(0)). (11) 


too 


[tama 
Finally, using (9) and (11), we have 


lim (s Fa(s)) — f(0) = lim (f(t) — f(0)). 


s—0 too 
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Neither term on the right-hand side depends on s, so we can remove the limit 
and simplify, resulting in the final value theorem 


3 Stability analysis of linear conformable fractional 
differential system 


In this section, we consider the stability of the following linear conformable 
fractional differential system 


tT, u(t) = Ax(t), t>0, x(0)=20, (12) 
where x € R®, matrix A € R"*", ro = (x10,220,---,2no)?, and a = 
[01, @2,...,Q»] such that 0 < a; <1, fori =1,2,...,n. 

Remark 1. If a = a] = ag =-:: = Qn, then system (12) is called a com- 


mensurate order system, otherwise system (12) indicates an in-commensurate 
order system. 


Definition 4. The zero solution of linear conformable fractional differen- 
tial system (12) is said to be stable if, for any initial value xo, there exists 
an € > 0 such that ||a(t)|| < € for all t > to. The zero solution is said to be 
asymptotically stable if, in addition to being stable, ||x(t)|| > 0 as t — oo. 

Now we state stability theorems from commensurate order system, next 
produce conditions for asymptotical stability of the in-commensurate order 
system. 


Theorem 6. The solution to the linear commensurate order system (12) 
is given by ; 
a(t) = ap e24*, (13) 


where the solution is assumed to be differentiable on (0,00). 
Proof. Taking the fractional Laplace transform of (12), by formula (7), we 


have that 
sX,.(s) — s%q = AXQ(s), 


so that 
X,(s) = xo(sI — A)". (14) 


Now, by applying the fractional Laplace transform of the fractional exponen- 
tial function, we obtain the claimed result. 
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Definition 5. If all the eigenvalues of A satisfy 


larg(X(A)| > 5. (15) 


then A is said to be Hurwitz stable. If all the eigenvalues of A satisfy 
larg(A(A)| 2 7/2, (16) 
and if |arg(Aj;0(A)| = 7/2, A;0 only correspond to simple elementary divisor 


of A, then A is said to be quasi-stable. 


Theorem 7.The zero solution of the conformable fractional differential sys- 
tem (12) is asymptotically stable if and only if A is a Hurwitz stable, and the 
zero solution of system (12) is stable but not asymptotically stable if and only 
if A is quasi-stable. 


Proof. According to the relation (13), we have 
x(t) = x ea At 


Let A = SJS—!, J is a Jordan canonical form, where J;, 1 < i <r has the 
following form 


ri 1 
Me 
eS oa ; rx EC, 
Ai 1 
Xi NXni 
and )>;_, nj =n. Then, 
als a 1 —lya a a 
eat = ea rds mw _g eavt s-4, 
1 yyo ; a ae a a @ 
eat = diag (eo seat eadrt Me 
tO oes i ake RR A ke a 
a a2 2! aia! (nj—1)! 
ise 1 pty Be 
: 0 1 ail ai? (nj—2)! ; 
ea lit = 00 . earst : 
; lia 
i. 1 at 
0 0 :-- O 1 abe 


One can easily show that the stability of zero solution of system (12) is 


determined by the boundedness of e«4**, or the boundedness of ez7*" or the 


20 H. Rezazadeh, H. Aminikhah, and A.H. Refahi Sheikhani 


boundedness of all eas!" j = 1,2,...,7r ie., Jarg(A;(A)| > 7/2, and when 
Jarg(A;(A)| = 7/2, nj; = 1, ie., A is quasi-stable. Asymptotical stability of 
the zero solution of system (12) is given by 


lim e=4# 0, 
t—+00 
or 
: Ly 
lim ex%* 0, 
t—+00 


which is equivalent to 


lim e#73**°=0, (j=1,2,...,r). 


t—+00 
That is, 
1 Li 47 jie ke ee 
1 ae a? 2. °° *° ars—t (nj—D! 
1 1 £6r 5 —2)0 
0 if tb a2 (n;—2)! 
: . 1) 4° : 
jim | 0 0 eh eee : eo it = 0, Wp S12 tj 
: : 1 
a an oe 


NZXN; 


and thus |arg(\;(A)| > 7/2 (j = 1,2,...,n) implying that A is a Hurwitz 
matrix. The proof is complete. 


Remark 2. If there is a Xo such that |arg(Ao)| < 7/2, then zero solution of 
system (12) is unstable. 


Remark 3. If A has zero eigenvalue, then zero solution of system (12) 
is unstable. 


Proof. If \ = 0, then from the proof of Theorem 7, we have 


{G-Da 


tI-Deo 


GEG = treme Up eae 


A=0 


It is obvious that jim, Gepaa=n = © for 7 > 1. Thus, jim, ||a(t)|| = co. 


Now, we consider an in-commensurate linear conformable fractional dif- 
ferential system 


tT, x1(t) == a1121(t) + a1222(t) ya ae aa Gin&n(t), 
tTovy L2(t) = ao101(t) + aeate(t) +--+: + dantn(t), 


(17) 


tTa,, Un(t) = Ani @1(t) + Qne%o(t) + +++ + Anntn(t), 
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where x;(0) = vj9 andO<a; <1 fori=1,2,...,n. 
We study the stability of system (17) by applying the fractional Laplace 
transforms on both sides of this system, we have 


8X q; (s) — Lio = x dij Xai (s), (18) 
j=l 
fori =1,...,n, where X,,(s) is the fractional Laplace transform of x;(t). 


We can rewrite (18) as follows 


Xa, (s) 
Xa (8) 
A(s) : = Xo. (19) 
Xa, (8) 
in which 
Aii(s) Aio(s) . Ain(s) 
Aoi (s) Ao0(s) a Aon(s) 
A(s = : ) 
Ani(s) Ano(s) Ann(s) 
where 
vy oe S— Ay ifi= J; 
Auj(s) = { ai; otherwise. 
and 29 = (10,%20,---;%no)". For simplicity, we call A(s) a characteristic 


matrix of (17), moreover det(A(s)) = 0 is the characteristic equation of sys- 
tem (17). Now, we express the main theorem for checking the stability of 
system (17). 


Theorem 8. /f all roots of det (A(s)) = 0 have negative real parts, then 
zero solution system of (17) is asymptotically stable. 


Proof. Multiplying s on both sides of (19) gives, we have 


A(s). = 8X0. (20) 
sXo,,(S) 


if all roots of the det (A(s)) = 0 lie in open left half complex plane (i.e. 
R(s) < 0), then, we consider (20) in R(s) > 0. In this restricted area, the 
relation (20) has a unique solution sX(s) = (sXq, (8), $Xqa.(S),.--,;5Xa, (s)). 
Since lim s = 0, so we have 
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lim  _sXq,(s) = 90, P= 125M 
80, R(s)>0 


which from the Theorem 5, we get 


lim a(t) = lim (a1(t), vo(t),...,@n(t)) = lim (sXq,(s), 8Xa,(8),...,$Xa, (s)) = 0. 


t—0o t-0co s—0 


The above result shows that system (17) is asymptotically stable. 


The inertia of a matrix is the triplet of the numbers of eigenvalues of A 
with positive, negative, and zero real parts. Now, we generalize the inertia 
concept for analyzing the stability of fractional linear system. 


Definition 6. The inertia of the system (17) is the triple 


In(a)(A) = (Beste) (A), Un(a) (A), On(a) (A)), 


where Tp(a)(A), Yn(a)(A) and dbn(a)(A) are, respectively, the number of roots 
of det (A(s)) = 0 with positive, negative, and zero real parts. 


Theorem 9.The linear conformable fractional differential system (17) is 
asymptotically stable if any of the following equivalent conditions holds. 

(i) The matriz A is a Hurwitz matric. 

(it) Tn(a) (A) = bn(a) (A) = 0. 

(itt) All roots of the characteristic equation of system (17) satisfy |arg(s)| > 
m/2. 


Proof. According to Theorem 8 and Definition 5, proof can be easily ob- 
tained. 


A Stability of non-linear conformable fractional 
differential systems 


In this section, we will mainly discuss the stability of a nonlinear conformable 
fractional differential system, which can be described by 


:T,, X(t) = F(X(t)), t>0, 0O<a<1, (21) 
with the initial value X(0) = Xo, where 


fi(ai(t), vo(t),...,¢n(t)) 


fo(ai(t), vo(t),..., @n(t)) 
rex = | RDO 


Pace Ons aR) 
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and X(t) = (x1(t), zo(t),...,2n(t))? € R”. 

Theorem 10. Let X(t) = (25 (0); 23 (0); 5235 28n(2))* is the equilibrium of 
system (21), ie. 4TaX = F(X) = 0 and J= (22) | px is the Jacobian 


matrix at the point X, then the point X is asymptotically stable if and only 
if J is a Hurwitz matric. 


Proof. Let C(t) = X(t) — X, where C(t) = (4:(t), Go(t),...,¢n(t)) is a small 
disturbance from a fixed point. Therefore 


Ta C(t) = Ta (X()-X), (22) 
since +Ty (xo _ x) = ,T,X(t) — 1T»X, and :T,X = 0, thus, we have 


tTo C(t) = To X(t) = F(X(t)) = F(t) + X) 


= F(X) + JC(t) + higher order terms 
x JC(t). 
System (22) can be written as 


tT C(t) © JC(t), (23) 


with the initial value ¢(0) = Xo — X. The analytical procedure of lineariza- 
tion is based on the fact that if the matrix J has no purely imaginary eigen- 
values, then the trajectories of the nonlinear system in the neighborhood 
of the equilibrium point have the same form as the trajectories of the lin- 
ear system [17]. Hence, by applying Theorem 7 the linear system (23) is 
asymptotically stable if and only if all roots of the characteristic function 
of J satisfy |arg(A(J)| > 7/2, which implies that the equilibrium X of the 
nonlinear conformable fractional system (21) is as asymptotically stable. 


Remark 4. The nonlinear conformable fractional system (22) in the point 
X is asymptotically stable if and only if ta) (J) = dna) (J) = 0. 


5 The multi-step homotopy perturbation method 


The homotopy perturbation method (HPM) was proposed by He [10] in 1999. 
This method has been used by many mathematicians and engineers to solve 
various functional equations. Although the HPM yields a solution series 
which converges very rapidly in most linear and nonlinear equations, in the 
case of a large time interval t it may produce a large error. To overcome this 
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shortcoming, the multi-step homotopy perturbation method (MHPM) was 
presented in [5], to solve nonlinear ordinary differential equations. For the 
convenience of the reader, we will first present a brief account of HPM. Let 
us consider the following differential equation: 


A(u) — f(r) =0, reQ, (24) 


with boundary conditions 
aa 0, rE Tr, (25) 


where A is a general differential operator, B is a boundary operator, f(r) is 
a known analytic function, and [ is the boundary of the domain 2. 
The operator A can be generally divided into two parts LZ and N, where L 
is linear, while N is nonlinear. Therefore Equation (24) can be written as 
follows: 

L(u) + N(u) — f(r) =0. (26) 


By using homotopy technique, one can construct a homotopy y(r,p) : 
Q x [0,1] —> R which satisfies 


A(y,p) = 1—p)[L(y)+ Luo) +PIA(y)—F(r)] = 9, pe [0,1], re Q, (27) 


which is equivalent to 


H(y,p) = L(y) — L(uo) + pL(uo) + pIN(y) — f(r)] = 0, (28) 


where p € [0,1] is an embedding parameter, and uo is an initial guess approx- 
imation of Equation (24) which satisfies the boundary conditions. Clearly, 
we have 


I 


HA(y,0) = L(y) — L(uo) 


0, (29) 
H(y,1) = Ay) — f(r)] = 0. 

Thus, the changing process of p from 0 to 1 is just that of y(r,p) from 
uo(r) to y(r). In topology this is called deformation and L(y) — L(ug) and 
A(y) — f(r) are called homotopic. If, the embedding parameter p, (0 < p< 
1) is considered as a small parameter, applying the classical perturbation 
technique, we can naturally assume that the solution of Equations (27) and 
(28) can be given as a power series in p, i.e., 


y=yotpy tp yet. (31) 


According to HPM, the approximation solution of Equation (24) can be ex- 
pressed as a series of the power of p, i.e. 


u=lmy=ytyityet--:. (32) 
pol 


Stability Analysis of Conformable Fractional Systems 25 


Now, consider a general systems of conformable fractional ordinary dif- 
ferential equations 


tT, U1 a gilt, U1, U2, tee Um) am fi(t), 


tT. U2 + go(t, U1, U2,...,Um) = fo(t), (33) 


Paes + Gm(t, U1, U2,---,Um) = fin (8), 
subject to the initial conditions 
ui(to) =c1, Ueal(to) = ce, .-. Um(to) = Cm. (34) 
First, we write the system (33) in the operator form 


L(u1) 2 Ni(t, U1, Ua, . .+)Um) _ fi(t) — 0, 


L(uz) + No(t, ui, u2,.--,;Um) — fo(t) = 0, (35) 


L(um) + Nm(t, U1, Ua, ate Um) = fim(t) J 0, 


subject to the initial conditions (34), where L (uw) =: Tq (uz) is a linear oper- 
ator and Ni, No,..., Nm are the nonlinear operators. To apply the MHPM, 
we first construct a homotopy for system (35) as follows 


L(ui) — L(v1) + pL(v1) + p[Ni (ui, ua,---,Um) — fil] = 0, 


L(uz) — L(v2) + pL(v2) + p[No(ur, u2,..-,Um) — fo(t)] = 0, (36) 


L (ttm) — L(¥m) + pL(Ym) + plNmn tia t2, - stm) — fm (2)] = 0, 


where V1, U2,.--,Um are initial approximations which satisfying the given con- 
ditions. Let us take the initial approximations as follows 


u1,0(t) = v1(t) = ui(to) = «1, 


Ug,0(t) = 2(t) = Ua(to) = ea, (37) 


and 
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ur(E) => u1,o(t) = pur,1(t) a a pu1.2(t) TT p>u1,3(t) eaey 


ur(E) = u2,0(t) =e pugz,1(t) am p?u2,2(t) cr ied p>u2,3(t) eaey (38) 


Um(t) = Umo(t) + pUm,1(t) + p?um,2(t) 7 p?'Um,3(t) ney 


where u;,;, (¢ =1,2,...,m;7 =1,2,...) are functions yet to be determined. 
Substituting (38) into (36) and arranging the coefficients of the same powers 
of p, we get 


L(ui,1) + L(v1) + Ni (ut,0, U2,0,---5Umo) — fi =0, — U1,1(to) = 9, 


L(u2,1) + L(v2) + No(u1,0, U2,0,-.-,Um,o) — f2 =0, u2,1(to) = 9, 


L(Um,1) + L(vm) + Nm (1,0; U2,05++5 ,Um,0) <7 tim = 0, Um,1(to) = 0, 


(uae) + Ni (ua tasty, Um) = 0, taa(to) = 0, 


D(ua,2) + Ny (u1,1, U215-+-5Um1) =0,  Ue,2(to) = 0, 


L(um,2) + Nm(u1,1, U2,15+++;Um1) =0,  Um,2(to) = 0, 
(39) 
etc. We solve the above systems of equations for the unknowns w;,;, (¢ = 
1,2,...,m;j =1,2,...) by applying the inverse operator 


Pe if (.) da(t). (40) 


Therefore, according to HPM the n-term approximations for the solutions of 
(33) can be expressed as 


din (t) = u(t) = limp 1 wi (t) = RZ 6 ure (O), 


gayn(t) = Uo(t) = limp +1 u2(t) = x56 u2,(C), (41) 


dm.n (¢) = un (t) = limp+1 Um (i) = ‘aa Um,k(t), 


The solution obtained by HPM is not valid for large t. A simple way of 
ensuring validity of the approximations for large t is to treat the algorithm 
of HPM in a sequence of intervals choosing the initial approximations as 
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uro(t) = v1(t) = ui(t*) = cj, 


ug,o(t) = v2(t) = ua(t*) = 3, (42) 


ee) a Um(t) = Um (t*) = Cs 


where t* is the left-end point of each subinterval. Now we solve (42) for the 
unknowns u;;, (i = 1,2,...,m;j = 1,2,...) by applying the inverse linear 
operator 


* 


LG i (.) da(t). (43) 


In order to carry out the iterations in every subinterval of equal length 
At, (0, t1), [t1, ta), 
[t2,t3),..., [ty-1,t), we would need to know the values of the following 


wz o(t) = u(t"), uo(t) = ua(t"),---,Um,o(t) = Um(E*). (44) 


But, in general, we do not have these information at our clearance except 
at the initial point t* = tp. A simple way for obtaining the necessary values 


could be by means of the previous n-term approximations ¢1,n, 2,n,---; m,n 
of the preceding subinterval given by (41), i.e. 
tio Pinlt"), Uoo & Gan(t"),---5Umo = Omn(t). (45) 


6 An illustrative example 


The following example is presented to illustrate the effectiveness and appli- 
cability of the proposed stability criteria. 

The Lotka-Volterra equations, also known as the predator-prey (or parasite- 
host) equations, are a pair of first order, non-linear, differential equations 
frequently used to describe the dynamics of biological systems in which two 
species interact on each other, one is a predator and the other is its prey. They 
were proposed independently by Alfred J. Lotka in 1925 and Vito Volterra in 
1926 [9,16]. The conformable fractional Lotka-Volterra system is described 
by the following nonlinear fractional conformable differential equations [13] 


tT, ti(t) = 21 (t) (r —a ai(t) — b xo(t)), 
0<ay,a2q < 1, (46) 
tTa.t2(t) = £2(t)(-—d+ezi(t)), 


with the initial values x ,(0) = x19 and x2(0) = 229 , where 21,22 > 0 are 
prey and predator densities, respectively, and all constants r, a, b, c and d 
are positive. This system has three equilibrium E, = (0,0), Bz = (£,0) and 
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E3 = (4, 5%). The Jacobian matrix for equilibria E* = (2*,y*) is defined 


c? c 


0.85 


0.75 


Figure 1: The equilibrium point E2 of the system (46) with (a1,a2) = (0.95, 0.9) and 
(r, a, b, c,d) = (1,3, 1,2, 2) is asymptotically stable 


oe" Aa by* ba* 


cy —d+ca* (47) 


When (r,a,b,c,d) = (1,3,1,2,2) the system equilibrium points are FE, = 
(0,0), By = (%,0) and E3 = (1,—2). Corresponding eigenvalues for equi- 
librium point £; are A; = 1, Az = —2, since In(a)(J) = (1, 1,0), hence the 
equilibrium point FE) is unstable. For equilibrium point F2 the eigenvalues are 
Ay = —1,A2 = =, because In(a)(J) = (0,2,0), thus the equilibrium point 
is asymptotically stable and the equilibrium point EF; is unstable, because 
In()(J) = (1, 1,0). Let us consider the following parameters of the system 
(46), (r, a, b, c,d) = (3, 0.5, 2,2, 2.5), for these parameters the system (46) has 
three equilibrium points FE, = (0,0), Eg = (6,0) and F3 = (1.25, 1.1875) 
and their corresponding inertias are In(q)(J) = (1,1,0) for Ey, E2 and 
In(o)(J) = (0,2,0) for £3, thus the equilibrium point 3 is asymptotically 
stable. 

Figure 1 shows that system (46) with parameters (r, a, b,c, d) = (1,3, 1, 2,2) 
and (a1, @2) = (0.95,0.9) is asymptotically stable in the equilibrium point 
E>. From Figure 2 we can see that the system (46) is asymptotically stable 
in the equilibrium point F2, with parameters (r, a, b,c,d) = (3, 0.5, 2,2, 2.5) 
and (a1, Q2) = (0.95, 0.85). 
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34 


25-7 


0.55 


Figure 2: The equilibrium point E3 of the system (46) with (a1, a2) = (0.95,0.85) and 
(r, a, b, c,d) = (3, 0.5, 2,2, 2.5) is asymptotically stable 


In system (46) when we take a = 0, we obtain modified the fractional 
conformable Lotka-Volterra system 


tT, 41 (t) = ravi (t) — b x1 (t)xo(t), 
(48) 
tTo,%2(t) = £2(t)(—d+cx,(t)). 
System (48) has two equilibrium points: F, = (0,0), BH = (¢, ae 
The Jacobian matrix of the system (48), evaluated at the equilibrium 
E* = (a*,y*), is given by 


jai" by" —bx* 


cy —d+cx* (49) 


The eigenvalues of the Jacobian matrix (48) evaluated at all equilibrium 
points show that all equilibria are unstable. Since, for the equilibrium point 


FE, we obtain A; = r and Ag = —d, for the equilibrium point E2 we get 
Ai,2 = £ivVrd. All these eigenvalues satisfy the condition for the system to 


be unstable (In(a)(J) = (1,1,0)). Figure 3 shows that the system (48) with 
parameters (r, b,c, d) = (1,1,4,2) and (a1, a2) = (0.97, 0.97) is unstable. 

All the results are calculated by using the computer algebra package 
Maple. The term-number of MHPM series solutions is fixed N = 4 and the 
time step size h = 0.1, with the initial conditions (#1(0),x2(0)) = (0.2, 0.8). 
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2.57 


0.55 


Figure 3: System (48) with (a1, a2) = (0.97, 0.97) and (r, b,c, d) = (1,1, 4, 2) 
is unstable 


7 Conclusion 


In this paper, we have studied the stability analysis of fractional differen- 
tial systems. Fractional derivatives are described by conformable fractional 
derivatives. At first, the stability conditions established by fractional ex- 
ponential function for commensurate linear fractional differential systems. 
Then, we proved a fractional version of final value theorem and by using this 
theorem we proposed sufficient conditions on the asymptotical stability for 
in-commensurate linear fractional differential system. The numerical simula- 
tions of conformable fractional Lotka-Volterra system are used to illustrate 
our main result. Although this paper just focuses on the linear systems and 
non-linear systems but, the problems remain open for the multi-order frac- 
tional differential systems, the time-delayed fractional differential systems 
and distributed order fractional systems. This will be the investigation goal 
of future works. 
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